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Abstract 


g 

Cfl 

o 

We use the NASA GEOS Chemistry Climate Model (GEOSCCM) to quantify the con- 
tribution of two most important brominated very short-lived substances (VSLS), bromo- 
form (CHBr 3 ) and dibromomethane (CH 2 Br 2 ), to stratospheric bromine and its sensitiv- 
ity to convection strength. Model simulations suggest that the most active transport of 
VSLS from the marine boundary layer through the tropopause occurs over the tropical 
Indian Ocean, the Western Pacific warm pool, and off the Pacific coast of Mexico. To- 
gether, convective lofting of CHBr 3 and CH 2 Br 2 and their degradation products supplies 
~ 8ppt total bromine to the base of the Tropical Tropopause Layer (TTL, ~ 150 hPa), 

( f) 

similar to the amount of VSLS organic bromine available in the marine boundary layer 
(~ 7.8-8.4ppt) in the above active convective lofting regions. Of the total ~ 8 ppt VSLS- 
originated bromine that enters the base of TTL at ~ 150hPa, half is in the form of 
source gas injection (SGI) and half as product gas injection (PGI). Only a small portion 
(< 1 0 %) the VSLS-originated bromine is removed via wet scavenging in the TTL before 
reaching the lower stratosphere. On global and annual average, CHBr 3 and CH 2 Br 2 , to- 
gether, contribute ~ 7.7pptv to the present-day inorganic bromine in the stratosphere. 
However, varying model deep convection strength between maximum and minimum 
convection conditions can introduce a ~ 2.6 pptv uncertainty in the contribution of VSLS 
to inorganic bromine in the stratosphere (Br y ). Contrary to the conventional wisdom, 

minimum convection condition leads to a larger Br y as the reduced scavenging in 
soluble product gases, thus a significant increase in PGI (2-3 ppt), greatly exceeds the 
relative minor decrease in SGI (a few lOths ppt). 

g 

o 

1 Introduction 

o' 

Very short-lived (VSL) bromocarbons originate mostly from ocean biogenic sources, 
and when transported into the stratosphere, they exert a significant impact on the 
bromine budget and stratospheric ozone depletion (Kurylo and Rodriguez, 1999; 
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Sturges et al., 2000). Recent years have seen significant progress in modeling ef- 
forts to quantify the contribution of brominated VSL substances (VSLS) to strato- 
spheric inorganic bromine (Br y ) (e.g. Warwick et al., 2006; Liang et al., 2010; 
Aschmann et al., 2011; Aschmann and Sinnhuber, 2013; Hossaini et al., 2012a, b, 
2013). These modeling studies suggest that brominated VSLS contributes 4.8-7 ppt 
to reactive stratospheric bromine, in good agreement with the 1-8 ppt estimate from 
satellite and balloon-borne observations (WMO 201 1 ; Sinnhuber et al., 2005; Sioris et 
al., 2006; Dorf et al., 2006a, 2008; Salawitch et al., 2010). 

The most important VSLS pathway to the stratosphere is via convective lofting 
through the tropical tropopause layer (TTL). Atmospheric chemistry and transport of 
brominated VSLS involves the coupling of various complex processes, e.g. highly un- 
uniform ocean emissions, convective transport of source gases and product gases, as 
well as wet scavenging of soluble product gases. An accurate modeling representation 
of these processes remains a challenging task. Modeling of these processes in general 
requires some extent of simplified assumptions and/or parameterizations, which may 
differ greatly from one model study to another. These differences have lead to signifi- 

VSLS 

cant differences in the modeled Br y estimates. For example, a recent model study 

\ IQI Q 

by Hossaini et al. (2013) found that modeled Br y could vary by a factor of 2 when 
using four different recently published emission inventories (Liang et al., 2010; Pyle et 
al., 2011; Ordonez et al., 2012; Ziska et al., 2013). In general, model studies agree 
relatively well on source gas injection (SGI), suggesting about ~50% of bromoform 
(CHBr 3 ) and ~ 90% of dibromomethane (CH 2 Br 2 ) can reach the stratosphere through 
SGI (e.g. Dvortsov et al., 1999; Nielsen and Douglass, 2001; Aschmann et al., 2011; 
Hossaini et al., 2010, 2012a; Brinckmann et al., 2012). However, modeled contribution 
on product gas injection (PGI) is highly uncertain depending on how wet scavenging 
is implemented. Many early model studies assume a uniform washout lifetime against 
wet scavenging (e.g. Dvortsov et al., 1999; Nielsen and Douglass, 2001; Hossaini et 
al., 2010; Aschmann et al., 2009), which was inadequate and lead to an underestimate 
in modeled Br y SLS (Hossaini et al., 2012a; Aschmann et al., 2011). A recent study by 
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WQIQ 

Aschmann and Sinnhuber (2013) shows that treating Br y as a single soluble tracer 

WQIQ 

(e.g. Liang et al., 2010) can also lead to an underestimate. They found that Br y 
increased from 3.4 ppt to 5ppt when switching from an idealized setup with a single 
soluble inorganic bromine tracer to a full chemistry scheme. Thus, an accurate mod- 
5 eling of transport and wet scavenging of PGI seems to be the key in narrowing the 

VSLS 

uncertainty of model estimate of Br y 

Understanding how the contribution of VSLS to stratospheric bromine varies with 
convection strength has significant climate implications. A recent modeling analysis 
by Hossaini et al. (2012b) suggests that VSLS source gas injection increases from 
io ~ 1.7 ppt in 2000 to ~ 2. 0-2. 7 ppt in 2100 using the Intergovernmental Panel on Cli- 
mate Change (IPCC) representative concentration pathways (RCPs) scenarios, as 
the future simulations feature stronger tropical deep convection transport to the lower 

WQIQ 

stratosphere. On the other hand, the overall response of Br y to convection strength 
is somewhat murky. Earlier studies deploying a uniform washout lifetime found differ- 
i 5 ent washout rates result a significant range in the contribution of CHBr 3 to Br y , 0.5- 
3 ppt in Sinnhuber and Folkins (2006) and 1.6-3 ppt in Aschmann et al. (2011). Liang 
et al. ( 2010 ) which deploys explicit wet scavenging in convective updrafts found that 
convective scavenging only accounts for ~ 0.2 ppt difference in modeled Br y . Sim- 
ilar results were reported in Aschmann et al. (2011), who showed while SGI is highly 
20 correlated with convective activity, the impact on total stratospheric bromine in a full 
chemistry scheme is nearly insensitive to dehydration, likely due to convection dilution 
and increased scavenging. 

In this study, we use a 3-dimensional Chemistry Climate Model with fully interactive 
CHBr 3 and CH 2 Br 2 , the two most important brominated VSLS, to better understand 
25 how VSLS and their degradation products enter stratosphere. We will also test the 
uncertainty in modeled Br y due to varying strength in convection and scavenging. 
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2 Model 


We conduct model simulations using the NASA GEOS Chemistry Climate Model 
(GEOSCCM) Version 2, which couples the GEOS-5 GCM (Reinecker et al., 2008) with 
a stratospheric chemistry module. Model simulations have a horizontal resolution of 
5 2.5° x 2° (longitude by latitude) with 72 vertical layers from surface to 0.01 hPa. The 
GEOSCCM model uses a flux-form semi-Lagrangian dynamical core (Lin, 2004) and 
the Relaxed Arakawa Schubert (RAS) parameterization for convection (Moorthi and 
Suarez, 1992). The GEOS-5 moist processes are represented using a convective pa- 
rameterization and prognostic cloud scheme. All chemical kinetics and photolysis rates 
io are calculated following JPL 2010 (Sander et al., 2011). 

The GEOSCCM V2 stratospheric chemistry module includes all important gas phase 
stratospheric reactions as described in Douglass and Kawa (1999). We modify the 
chemistry scheme to include CHBr 3 and CH 2 Br 2 . The two organic source gases are re- 
leased at the ocean surface following the geographically resolved emission distribution 
is described in Liang et al. (2010) and are destroyed in the atmosphere via photolysis and 
reaction with the hydroxyl radical (OH). OH above the tropopause is calculated online 
in the stratospheric chemistry module. OH in the troposphere is relaxed to a monthly 
mean climatological fields documented in Spivakovsky et al. (2000). Global annual 
mean OH from Spivakovsky et al. (2000) is 1.16 x 10 6 molecules cnrT 3 , yielding an at- 
20 mospheric methyl chloroform (CH 3 CCI 3 ) lifetime of ~5.5yr. Below tropopause, only 
simple bromine chemistry is considered which partition 80 % of inorganic bromine prod- 
ucts into HBr and the remaining 20% as HOBr, ratios adopted from Yang et al. (2005). 
Above the tropopause, bromine from degradation of CHBr 3 and CH 2 Br 2 is released 
as Br and interacts fully with stratospheric chemistry. We assume high solubility for 
25 all three inorganic bromine reservoir species, HBr, HOBr, and BrON0 2 , and compute 
their wet scavenging. The wet removal processes are detailed in Liang et al. (2010), 
which include scavenging in rainout (in-cloud precipitation) and washout (below-cloud 
precipitation) in both large-scale precipitation (Giorgi and Chameides, 1986) and deep 
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convective updrafts (Balkanski et al., 1993). We also consider re-evaporation back to 
the atmosphere. No heterogeneous chemistry for VSLS is included in our simulations, 
but results from Aschmann and Sinnhuber (2013) suggest that the impact of heteroge- 
neous chemistry is minor and the inclusion of heterogeneous activation prevents loss 
5 in scavenging and can increase Br y by 10%. 

In this study, we present results from two 51 yr simulations between 1960-2010, one 
with brominated VSLS chemistry (/? V sls) anc * one without (/? B ase)> to examine tropo- 
sphere to stratosphere transport (TST) of VSLS and their contribution to stratospheric 
bromine. It takes a considerable long time (~ 1 5 yr) for the full impact of VSLS to reach 
io the upper stratosphere, thus the first 1 6 yr are considered as spin-up. A full model eval- 
uation and discussion of the impact of VSLS on stratospheric ozone are presented in 
a separate paper, currently in preparation. For this study, we focus on model results 
from the last 1 0 yr of the two simulations. 

To quantify the uncertainty in Br y due to changes in convection strength, we con- 
15 duct two sensitivity simulations from 1980 to 2010 to represent minimum convection 
(R mincnv) ar *d maximum convection (#maxcnv) conditions by varying five convective 
parameters in RAS that impact the strength of deep convection, clouds, convective con- 
densate and re-evaporation (Ott et al., 2009, 2011). These five convective parameters 
are identified as the most strongly influencing out of the total 16 parameters examined 
20 using a large number of “Monte Carlo” type simulations and ensemble simulations in 
both a Single-Column Model as well as the GEOS-5 GCM (Ott et al., 2011). We aim 
to use these two extreme conditions to the range of uncertainty that can be introduced 
due to variations in the strength of convection and wet scavenging. 
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3 Results 

3.1 The contribution of VSLS to stratospheric bromine 

The simulated CHBr 3 and CH 2 Br 2 have been evaluated extensively in Liang 
et al. ( 2010 ) and compares well with aircraft and surface observations in their atmo- 
5 spheric distribution, geographically and vertically, as well as the associated seasonality. 
Here we present a summary comparison of observed and modeled CHBr 3 in the tro- 
posphere for 30-60° S, 30°S-30°N, 30-60° N, and 60-90° N latitude bands (Fig. 1). 
The observed profiles are compiled using Whole Air Sampler canister measurements 
(Schauffler et al., 1999; Blake et al., 2003) from eight NASA aircraft missions, as de- 
io tailed in Liang et al. (2010), and averaged within each latitude band at 1 km vertical 
interval. The comparison clearly shows that the model captures well the source gas 
concentrations and vertical gradients in the tropics, mid and high latitudes in both hemi- 
spheres. 

Figure 2 shows a comparison of our simulated BrO with balloon measurements from 
is the LPMA (Limb Profile Monitor of the Atmosphere )/DOAS (Differential Optical Absorp- 
tion Spectroscopy) at Teresina (5.1° S), Aire sur I’Adour (43.7° N), and Kiruna (67.9° N) 
collected between 2003-2005 (Dorf et al., 2006a, b, 2008; Rozanov et al., 2011). The 
lower levels of BrO from /? BASE ' n comparison to the measurements indicate the es- 
sential role of VSLS in completing the stratospheric bromine budget. With VSLS, the 
20 GEOSCCM model simulates well the observed BrO at all sampled locations in the trop- 
ics, mid and high latitudes. Since the model is running freely by specifying only surface 
source gas emissions, the fact that the model simulates well the observed concentra- 
tions of both the source gases and BrO suggests that the model presents a credible 
representation of stratospheric bromine chemistry. 

25 Comparing results from the /? VSLS an d ^base runs (Fig- 3), we find that the in- 
clusion of CHBr 3 and CH 2 Br 2 adds a uniform ~7.7ppt Br y throughout most of the 

WQIQ 

stratosphere. The current Br y estimate is 55% higher than our previous estimate of 
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~ 5 ppt (Liang et al., 2010) and ~ 6 ppt estimate in Hossaini et al. (2013), though the 
latter two are driven with the same VSLS emissions. Compared to the idealized case in 
Liang et al. (2010) that tracks Bry SLS in a single highly soluble tracer, Bry SLS is present 
in both soluble forms (HBr, HOBr, and BrON0 2 ) and insoluble forms (Br, BrO, BrCI) in 
5 the fully interactive chemistry scheme. The ratio of insoluble and soluble Br y varies with 
time of the day as well as altitude and location. At daytime, when most convective lofting 
occurs, the majority of Br y exists as insoluble Br and BrO, therefore greatly increases 

the amount of Br y that survives wet scavenging during TST and ultimately reaches 
the stratosphere. The increase from ~ 5 ppt to ~ 7.7 ppt when switching from an ideal- 
io ized case to fully interactive chemistry is consistent with results from Aschmann and 
Sinnhuber (201 3), who found that Br y increased from 3.4 ppt to 5 ppt when switching 
from an idealized setup with a single soluble inorganic bromine tracer to a full chem- 
istry scheme. The different estimates between Hossaini et al. (2013) and this work is 
likely due to how the ratio of soluble and insoluble inorganic bromine is determined 
is inside the two models. Compared to this work which partitions inorganic bromine in its 
various forms based on chemical reaction rates, Hossaini et al. (2013) used a mean 

VSLS 

altitude-dependent HBr : Br y ratio. This likely leads to an excessive washout of Br y 
as during daytime when most of the convective lofting occurs, the majority of the inor- 
ganic bromine exists as insoluble Br and BrO. This suggests that a close approximate 
20 of soluble and insoluble inorganic product gases and the associated diurnal variation 
are critical to accurately quantify the contribution of VSLS to stratospheric halogen in 
VSLS modeling. 

3.2 Troposphere-to-stratosphere transport of VSLS 

We use the simulated monthly mean CHBr 3 distribution on the 355 K potential tem- 
25 perature surface (just below TTL) to show important tropical regions where active TST 
initiates (Fig. 4). Theoretically, trace gas distribution on the 365 K surface should be 
more indicative of TST as 365 K marks the zero radiative heating and air mass ele- 
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vated above this level can enter the lower stratosphere through slow radiative ascend 
(e.g. Gettelman and Forster, 2002; Fueglistaler et al., 2009). However, it is difficult 
to identify active TST regions on the 365 K map (not shown) as a significant portion 
of CHBr 3 is converted to product gases. This is not surprising for a short-lived com- 
5 pound with lifetime ~ 26 days (WMO, 2011) while on average it takes about 10 days 
for air to transport by ±10 K (Fueglistaler et al., 2004; Levine et al., 2007). Compared 
to CHBr 3 , it is much more difficult to identify active TST regions on a CH 2 Br 2 map, 
as CH 2 Br 2 distribution appears more zonal with a smaller meridional gradient due to 
a longer lifetime (~ 120 days, WMO, 2011) and thus more mixing with the surrounding 
io background air (not shown). The 355 K CHBr 3 map shows three active regions that 
can efficiently deliver VSLS to the base of TTL: (1) the tropical Indian ocean, (2) the 
Western Pacific warm pool, and (3) off the Pacific coast of Mexico (Fig. 4). There is 
significant seasonality associated with each entry region. The Indian Ocean appears 
as the most active region for TST of CHBr 3 and occurs all year long with a maximum 
is in boreal winter (DJF). Lofting in the Western Pacific reaches its maximum in boreal 
summer (JJA), while TST off the Mexico coast occurs mostly in boreal summer (JJA) 
and fall (SON). The importance of convective lofting in the Western Pacific warm pool 
in TST has been noted in many previous studies (e.g. Hatsushika and Yamazaki, 2003; 
Fueglistaler et al., 2004; Aschmann et al., 2009; Hossaini et al., 2012a; Aschfold et al., 
20 2012) with several others suggesting that the Indian Ocean is also an important region 
in TST of VSLS (Levine et al., 2007, 2008; Brioude et al., 2010; Hoyle et al., 2011). 

We show in Fig. 5 the vertical profiles of organic source gases (CH 2 Br 2 x 2, CHBr 3 x 3) 
and inorganic product gases (Bry SLS = Br + BrO + BrCI + HOBr + HBr + BrON0 2 ) in the 
three critical convective lofting regions to illustrate the transport and washout of the 
25 brominated VSLS during TST and the relative importance of SGI and PGI. Over the 
tropical Indian Ocean, on annual average, the mixing ratio of total bromine from VSLS 
(organic + inorganic) at 150hPa (~ 355 K) is ~ 8.5 ppt, the same as its surface abun- 
dance (~ 8.5 ppt). This implies convective lofting in this region is so efficient that the 
amount of VSL bromine that enters the base of TTL is set by the concentration at the 
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WQIQ 

marine boundary layer. As air ascends to higher altitudes, a small portion of Br y 
(~ 0.8 ppt, 10 %) is gradually removed from the atmosphere via wet scavenging until 
reaching a constant 7.7 ppt above lOhPa. The somewhat wavy structure in the total 
bromine profile in the lower stratosphere implies that not all wet scavenging happen 
5 locally and the dip in total bromine concentration between 100-50 hPa is likely associ- 
ated with recirculation of air back from the mid-latitude lower stratosphere where addi- 
tional scavenging occurs via in-cloud rainout. Results from the Western Pacific and off 
the coast of Mexico are similar, with a smaller surface abundance of organic bromine 
(7.8 ~ 8.1 ppt) and slightly less washout around lOOhPa. The global averaged total 
io bromine from VSLS (Br ) profile shows a maximum value of ~ 8 ppt at ~ lOOhPa. 
To supply this amount of Br to the LS, transport from the MBL has to initiate from 
tropical regions where active convective lofting is co-located with high surface concen- 
trations, where the collective CHBr 3 and CH 2 Br 2 abundance exceeds 8 ppt Br. 

Mapping streamlines on top of the VSLS organic bromine distribution in the tropics 
is clearly illustrates the importance of co-location of deep convection with high surface 
concentration regions (Fig. 6 ). Among the three ascending branches of the Walker 
circulation (Webster, 1983), while ascending in the Western Pacific penetrates deep- 
est into the TTL, ascending in the Indian Ocean is capable in delivering more VSLS 
bromine to the base of TTL due to higher surface concentrations. The surface abun- 
20 dance of VSLS bromine is largely dependent on the emission distribution used. Note 
Liang et al. (2010) assume uniform zonal emission strength across all longitudes with 
a few surface observation constraints from the Western Pacific and none from the trop- 
ical Indian Ocean. Hossaini et al. (2013) conducted a recent model study comparing 
four independent brominated VSLS emission estimates. Of all four emission inventories 
25 examined (Liang et al., 2010; Pyle et al., 201 1 ; Ordonez et al., 201 2; Ziska et al., 2013), 
modeled CHBr 3 and CH 2 Br 2 concentrations using the emission estimates from Liang 
et al. ( 2010 ) compare very well against tropical observations from multiple years of 
National Oceanic and Atmospheric Administration/Earth System Research Laboratory 
(NOAA/ESRL) surface flask measurements and the HIAPER Pole-to-Pole Observa- 
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tions (HIPPO) measurements between 2009-2011 supported by the National Science 
Foundation (NSF). When compared against the aircraft measurements collected in the 
Western Pacific during the Stratospheric Ozone: Halogen Impacts in a Varying At- 
mosphere (SHIVA) campaign, Liang et al. (2010) yields the best CH 2 Br 2 (mean bias 
5 of + 0.2 ppt compared with observations), compared to the other three emissions, but 
the simulated CHBr 3 is ~ +0.88 ppt too high on average. Note the high bias in modeled 
CHBr 3 in the Western Pacific will likely lead to a high bias in the model Br y , although 
the lofting in Indian Ocean seems to be more critical in determining the total amount 
of VSLS bromine in the lower stratosphere. The model simulated high concentration 
io of brominated VSLS over the Indian Ocean is yet to be evaluated when more surface 
observations become available and the importance of the Indian Ocean in delivering 
higher amounts of VSLS bromine into TTL needs to be assessed. The simulated high 
surface concentrations of CHBr 3 and CH 2 Br 2 near tropical Central America were vali- 
dated with the NASA TC 4 and INTEX-B measurements (Liang et al., 201 0). This region, 
is compared to the above two, is much less efficient in delivering VSLS bromine into the 
TTL due to the relatively weaker ascending motion. However, this ascending branch 
can be important in particular seasons and in individual years as the Walker circulation 
moves in the east-west direction between different phases of El Nino Southern Oscil- 
lation (ENSO). For example, using trajectory calculations, Levine et al. (2008) found 
20 a clear shift in the TTL origin of air parcels from the Western Pacific and Indonesia to 
those from the Eastern Pacific and South America in El Nino years. 

Of the total VSLS bromine that enter the base of TTL, about half (~ 4 ppt) is in the 
form of inorganic product gases (Fig. 5) despite our simple assumption that all VSLS 
degradation products are released in the troposphere as highly soluble HBr and HOBr. 
25 A commonly adopted approach in many previous modeling studies is to prescribe the 
VSLS concentration at the base of TTL with observed organic source gas mixing ra- 
tios in the upper troposphere and tracks the subsequent chemistry and transport (e.g. 
Sinnhuber and Folkins, 2006; Aschmann et al., 201 1 ; Aschmann and Sinnhuber, 2013). 

VSLS 

Our result implies that Br y estimate from such approach is not complete as it misses 
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an important component associated with PGI, which is as large as SGI at the base of 
TTL in our model simulation. A second implication of this result is that if we were to 
use measurements from the TTL region to quantify the contribution of VSLS to strato- 
spheric bromine, it is necessary to make measurements of both organic and inorganic 
5 forms to fully account for the impact of VSLS on the atmospheric bromine budget. 

3.3 The impact of convection strength 

We examine the difference between the two convection sensitivity simulations, #mincnv 
and /?maxrun> t0 illustrate the impact of deep convection on Br y (Fig. 7). Differ- 
ence between these two simulations (/? M incnv ~^maxrun) suggests that in regions that 
io mattered most (the three tropical convection centers), surprisingly, the minimum con- 
vection condition is more favorable for TST of VSLS. While weaker convection slightly 
decreases SGI (a few tenths ppt), PGI increases significantly (2-3 ppt) due to less 
scavenging under minimum convection condition. TST of the longer-lived CH 2 Br 2 is 
less sensitive to convection strength compared to that of CHBr 3 because of its longer 
is lifetime. The significant increase in CHBr 3 and CH 2 Br 2 in the mid-latitude bands under 
minimum convection mainly reflects the compensating responses in large-scale de- 
scent in a general circulation model as a result of changing convection strength (Ott et 
al., 2011). The decrease of the descent does not impact TST of VSLS. 

Globally, the differences in convection strength and wet scavenging introduce a 
20 ~ 2.6 ppt uncertainty in Br y (~ 6.5 ppt in /? MAXRUN and ~ 9.1 PP* in in run) (Fig. 3). 
Although these two simulations represent extreme parameterization conditions, BrO 
from both simulations still falls within the ±2.5 ppt uncertainty range of the DOAS bal- 
loon observations (Fig. 2, Dorf et al., 2006b). 
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4 Conclusions 

We use the NASA Goddard 3-D Chemistry Climate Model, GEOSCCM, to quantify the 
contribution of brominated very short-lived substances (VSLS) to reactive stratospheric 
bromine and to test its sensitivity to the strength of deep convection and the associated 
5 scavenging. 

The inclusion of CHBr 3 and CH 2 Br 2 in a fully interactive stratospheric chemistry 
module adds ~ 7.7 pptv to present-day stratospheric inorganic bromine. The most ac- 
tive transport of VSLS from the marine boundary layer (MBL) through the tropical 
tropopause layer (TTL) occurs where high surface concentrations of VSLS co-locate 
io with deep convection centers: (1) the tropical Indian Ocean, (2) the Western Pacific 
warm pool, and (3) off the Pacific coast of Mexico. On annual average, almost all VSLS- 
originated bromine available in the MBL (~ 7. 8 - 8 . 4 ppt) at these active convective loft- 
ing regions enters the TTL, half in the form of source gas injection (SGI) and half in the 
form of product gas injection (PGI). After lofting above 150hPa (~ 355 K), the majority 
is of the VSLS-originated bromine survives TST and reaches the lower stratosphere, with 
only a small portion (~ 0.8 ppt, 10%) removed by wet scavenging. Our model results 
point to a clear need for more surface measurements in the active VSLS lofting regions, 
in particular the tropical Indian Ocean and the Western Pacific, for improved emission 
estimate as well as a better quantification of how much VSLS bromine is available 
20 before ultimately enters the stratosphere through TTL. 

Our current estimate of the contribution of CHBr 3 and CH 2 Br 2 to stratospheric 
bromine (Br y ) is higher than previous modeling estimates, mainly reflecting the 
differences in how transport and wet scavenging of product gases are treated in dif- 
ferent models (Hossaini et al., 2012a; Aschmann and Sinnhuber, 2013). Interestingly, 
25 even though the models are driven with same emission estimates (Liang et al., 2010; 
Hossaini et al., 2013) and/or with similar source gas abundance at the base of TTL 
(Aschmann et al., 201 1 ; Aschmann and Sinnhuber, 2013), this study still yields higher 
Br y for two main reasons. Modeling VSLS and their degradation products in full 
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chemistry scheme leads to a better representation of the partition between soluble 
and insoluble product gases and the associated diurnal variation, both are critical in 

WQIQ 

an accurate Br y estimate. Secondly, tracking PGI from the free troposphere to the 
base of TTL, which is as large as SGI at the base of TTL in our model simulation, is 
5 also important. This was not included in many earlier model studies (e.g. Aschmann 
et al., 2011; Aschmann and Sinnhuber, 2013), thus implying a missing component in 
these earlier model VSLS contribution estimates. 

Differences in model deep convection strength can introduce ~30% (6.6-9.2ppt) 
uncertainty in simulated Br y within the same convection scheme. It is important to 
io point out that this uncertainty is likely larger among different models or a same model 
but different versions when different convective parameterizations are used. Contrary 
to conventional wisdom that the VSLS impact on stratospheric bromine is larger un- 
der more intense deep convection, our simulations suggest that minimum convection 
condition is favorable for TST of VSLS due to reduced scavenging in soluble product 
is gases under weaker convection condition. The impact of convection strength on PGI 
greatly outweighs the impact on SGI with the change in PGI ~ ten times larger than 
that of SGI. 
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Fig. 1. Comparison between the observed (red lines with horizontal bars indicating one stan- 
dard deviation) and simulated (black lines with gray shading showing one sigma variability) ver- 
tical profiles of CHBr 3 in the troposphere. Observations are compiled using Whole Air Sampler 
canister measurements from eight NASA aircraft missions are averaged for 30-60° S, 30° S- 
30° N, 30-60° N, and 60-90° N latitude bands at 1 km vertical interval. Model is sampled at the 
same location as the observations in the corresponding month in 2010. 
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Fig. 2. Comparison of simulated BrO from the GEOSCCM run /? VSLS (black solid line with 1- 
sigma variance in gray shading) with balloon measurements (red line with 1 -sigma uncertainty 
in horizontal bars) from the LPMA/DOAS Spectrometers in the stratosphere. The black dash- 
dotted line shows simulated BrO from the run without VSLS (Abase)- Simulated BrO from the 
two convective sensitivity simulations are also shown (dotted line for maximum convection con- 
dition and dashed line for minimum convection condition). The model is sampled at the same 
month and latitude as the balloon measurements but only daytime profiles from all available 
longitude grid points are collected to calculate the la variance of simulated BrO. 
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Fig. 3. Model simulated global and annual mean BrO and Br y profiles from the simulation with 
VSLS (/?vsls< solicl |in es) and the base simulation without VSLS (/? B ase> dash-dotted lines) for 
Year 2010. The gray shadings indicate the spread of annual mean Br y in the /? VSLS run with 
the hatched areas indicate the minimum to maximum range and the filled shadings indicate la 
variance. Br y from the two convective sensitivity simulations are also shown (blue dashed and 
dotted line for minimum and maximum convection conditions, respectively). 
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Fig. 4. lOyr averaged (2001-2010) seasonal mean distribution of simulated CHBr 3 (unit ppt 
Br) at the 355 K potential temperature layer (just below tropopause) in DJF, MAM, JJA and 
SON. The black boxes outline the three regions where most active troposphere-to-stratosphere 
transport occurs. 
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Fig. 5. The contribution of organic bromine (CH 2 Br 2 x 2, CHBr 3 x 3) and inorganic bromine 

WQI Q 

(Br y ) from CHBr 3 and CH 2 Br 2 degradation to atmospheric bromine in three active TST re- 
gions (a) the tropical Indian Ocean (left), (b) the Western Pacific warm pool (middle), and (c) 

WQI O 

the Pacific Coast of Mexico (right). Global mean Br y profile is also shown (thick black dotted 
line). The model results are 1 0 yr annual averages from 2001-2010. 
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Fig. 6. Longitude-height cross-section of modeled total VSL organic bromine (CHBr 3 x 
3+ CH 2 Br 2 x 2) (color contours) and streamlines (black lines) in the deep tropics. The organic 
bromine and wind streams are 1 0 yr annual averages (2001-2010) between 10°S-10°N. Sur- 
face organic bromine abundance is shown as the black solid line at the bottom panel. 
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Fig. 7. The simulated annual mean difference between the /? M incnv and ^maxrun runs in or ' 
ganic bromine (3x CHBr 3 , 2x CH 2 Br 2 ), inorganic bromine and total bromine from VSLS be- 
tween 355-380 K. The model results are lOyr annual averages from 2001-2010. 
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